Finite Temperature Collapse of a Bose Gas with Attractive Interactions 



0^ 



o 
O 

(N 



I 

o 
o 



> 
m 
m 

oo 

o 



X3 

o 
o 



X 




Normal 



100 10^ 10"^ 
Temperature [nK] 



10' 



Erich J. Muelleifl and Gordon Baymi 
Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urhana, IL 61801 

(October 26, 1999) 

We study the mechanical stabihty of the weakly interacting Bose gas with attractive interactions, 
and construct a unified picture of the collapse valid from the low temperature condensed region to 
the high temperature classical regime. As we show, the non-condensed particles play a crucial role 
in determining the region of stability. We extend our results to describe domain formation in spinor 
condensates. 
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The trapped Bose gas with attractive interactions is 
a novel physical system. At high densities these clouds 
are unstable against collapse; however at low densities 
they can be stabilized by quantum mechanical and en- 
tropic effects. Such stabihty and the subsequent collapse 
has been observed in clouds of degenerate ^Li Q. The 
collapse of the Bose gas shares many features of the grav- 
itational collapse of cold interstellar hydrogen, described 
by the Jeans instability related instabilities occur in 
supercooled vapors. Theoretical studies of the attrac- 
tive Bose gas, typically numerical, have been limited to 
zero or very low Q temperature. Here we give a sim- 
ple analytic description of the region of stability and the 
threshold for collapse, valid from zero temperature to well 
above the Bose condensation transition, and thus provide 
a consistent global picture of the instability. As we show, 
the non-condensed particles play a signficant role in the 
collapse at finite temperature Indeed in the normal 
state, collapse occurs even when no condensate exists. 

Our results are summarized in the phase diagram in 
Fig. |l|, which shows three regions: normal, Bose con- 
densed, and collapsed. This third region is not readily 
accessible experimentally; the system becomes unstable 
at the boundary (the solid line in the figure). The attrac- 
tive interactions, which drive the instability, have a char- 
acteristic energy V — (asn) / m per particle, where h is 
Planck's constant, {— — 1.45nm for ^Li) is the s-wave 
scattering length, n is the density, and m is the particle 
mass. At low temperature the only stabilizing force is 
quantum pressure, arising from the zero-point motion of 
the particles, with a characteristic energy Eq = ?i^/mi^, 
where L is the size of the cloud. At higher tempera- 
tures, thermal pressure, resulting from the decrease in 
entropy with decreasing volume of a system, also helps 
to stabilize the cloud. The characteristic energy scale for 
thermal pressure is Et = ksT, where T is the temper- 
ature. Naively, the line of collapse is where V becomes 
larger than both Et and Eq . Quantum pressure rapidly 
decreases as the system size increases, and in the thermo- 
dynamic limit the entire condensed phase is swallowed up 
by the collapsed phase (cf. Fig. 2). Due to the finite sys- 
tem size the line separating the normal and condensed 



FIG. 1. Phase diagram, density n [cm~'^] versus tempera- 
ture T [nK], for a small cloud of attractive bosons with spa- 
tial extent L — 3.15/im, scattering length Os = — 1.45nm, and 
atomic weight 7 (see [1]). Note the logarithmic scales. To ap- 
ply these results to harmonically trapped gases, the density 
n should be interpreted as the central density in the trap. 
The solid line separates the unstable (shaded) region from 
the stable region. 



phase is not a phase transition, but rather is a rapid 
cross-over. 

The zero temperature stability of a finite boson cloud 
can be understood in terms of the Bogoliubov excitation 
spectrum of a uniform gas [^: 



2 I k^V fc' 
w = I — +gn—, 

2m J m 



(1) 



where g = inh^as/m. In the attractive case, g < 0, all 
long wavelength modes with k"^ /2m < 2\g\n have imagi- 
nary frequencies and are unstable. A system of finite size 
L has modes only for k > 2tt/L and will have an exci- 
tation spectrum similar to Eq. (|l|) for short wavelengths. 
Thus one expects that a zero temperature attractive Bose 
gas will be stable if \g\n < h^n'^ /mL"^ We extend this 
analysis to finite temperature by calculating the zero fre- 
quency density-density response function xik). The on- 
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set of an instability in the longest wavelength mode is 
signaled by a divergence in = 27r/L). 

The main role of finite temperature, T, is to supple- 
ment quantum pressure with thermal pressure. In the 
limit T 3> Tc, quantum effects are negligible, and the line 
of collapse becomes the spinodal line of a classical liquid- 
gas phase transition, as characterized by Mermin To 
illustrate the physics of this non-degenerate regime we 
look for an instability in the uniform gas at zero wavevec- 
tor, fc = 0, calculating the compressibility in the Hartree 
approximation and effectively neglecting finite size 
effects such as quantum pressure. Within this approxi- 
mation the density is given by the self-consistent solution 
of 
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with Hartree quasiparticle energies e{p) — /2m + gn. 
Here fi is the chemical potential, and f3 — l/fc^T. In the 
classical limit, this equation becomes n = e^^'^~^"'> / Aj,, 
where the thermal wavelength is At = /mkBT)^/^ . 

The compressibility is proportional to the zero wavevec- 
tor density-density response function, x{k = 0) = 
—dn/dfi. We introduce the "bare" response, Xo{k = 0) = 
— (dn/dfj,)^, where the e are held fixed; in the classical 
limit Xo(0) = — /3n. The full response x(0), as in the 
random phase approximation (RPA) , is 



y(k) - ■ 

1 - 9Xa{k) 



(3) 



in the classical limit x(0) — + gPn). The col- 

lapse of the attractive system [g < 0) is signaled by a 
divergence of the compressibility. In the classical limit, 
this divergence occurs when \g\n — ksT [Q. The ap- 
proximations made here break down when ~ Yi' jml? 
and if applied outside their range of validity would er- 
roneously predict that the instability towards collapse 
prevents Bose condensation from occurring. 

In the experiments on ""Li, the atoms are held in 
a magnetic trap with a harmonic confining potential 
V{r) — \mijj^r^ ^ with w « 27r x 145s^"'^, and have a 
highly inhomogeneous density profile [Q. In connect- 
ing the theory of the uniform gas with the experiment, 
n should be interpreted as the density at the center of 
the trap, and the system size L as effectively ^Jh/mu) 
(= 3.2/im). As long as fcsT is larger than the trap en- 
ergy Tiuj « 7nK , the geometry of the trap should play no 
significant role in the collapse. 

We describe the system more generally using the re- 
sponse function derived from the shielded potential 
approximation pT| ], as used by Szepfalusy and Kondor 
to study the critical behavior of collective modes [ p^ . 
This approach, a simple example of the dielectric formal- 
ism, generates an excitation spectrum which is conserv- 
ing O] and gapless Q. At zero temperature it yields 



the Bogoliubov spectrum shown in Eq. (|^), and above 
Tc it becomes the standard RPA. Analysis in terms of 
Bogoliubov quasiparticles is left to a future paper. 

The susceptibility has the general form, Eq. (||), with 
the appropriate polarization part xo- Within the shielded 
potential approximation xa is the response of the non- 
interacting system. This response can be broken up into 
two parts, Xo ^-nd Xo; corresponding to the condensate 
and non-condensate contributions p5| , which at finite 
frequency are given by 
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Here no is the condensate density, Sk — k"^ /2m is the free 
particle spectrum, and the Bosc factors f{k) are given by 

(•g-/3(efc-Ai) _ 

Expanding Xo in the small parameter kAx we derive 
for T > Tc, 



gx'oik) = -2 
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where gv{z) = "^^z^ /j^ is the polylogarithm function. 
For chemical potential /i much larger in magnitude than 
/cbT, the system is in the classical regime, and Eq. ( |5a| ) 
reduces to gxo — "Pgn, as was found in the Hartree 
approach. Below Tc the chemical potential of the non- 
interacting system vanishes and the response functions 
are simpler: 



9Xl{k) = 
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n = no + ^C(3/2), 



(6a) 
(6b) 
(6c) 



with ({ly) — giy(l), the zeta function. Using these expres- 
sions we calculate the spinodal line separating the stable 
and unstable regions of Fig. ^ by setting k — 2t:/L and 
solving the equation 



X ' (X I - .g(xo + X'o) = 0. 



(7) 



Note that the condensate response Xo scales as L^, 
whereas below Tc the noncondensate response Xo scales 
as L. For realistic parameters, L is the largest length in 
the problem, so the condensate dominates the instability 
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FIG. 2. Scaling of the instability threshold with system 
size. Solid lines show the maximum stable density rimax for a 
given temperature. From the top, the system size is L = 0.3, 
3, 30 fim. The other parameters are the same as Fig. |^. The 
dashed line indicates the Bose-Einstein condensation transi- 
tion. Note the three scaling regimes; at low temperature, 



L , near the critical temperature Un 



L^/^ and 



at high temperatures n„ 
and (6)). 



is independent of L (see Eqs. (5) 



except when uq is much smaller than n. In Fig. |^ we plot 
the line of instability for various L. 

From the above equations we calculate the maximum 
stable value of the condensate density uq: 



4L3 



L^mkBT 



(8) 



('^■o)max decreases linearly with temperature, eventually 
vanishing at T = h'^/mL\a\. Conversely, the maximum 
total density increases rather quickly with temperature. 
Using parameters from the Rice experiments |l[ , we find 



i'K)max = 1706 -r/(8.8nif). 



(9) 



This result, plotted in Fig. |^, is consistent with the ex- 
periments, and agrees qualitatively with numerical mean- 
field calculations 0|. Quantitative agreement requires a 
more sophisticated treatment of the geometry. To verify 
the structure in Eq. and to map out the phase dia- 
gram in Fig. 0, future condensate experiments will need 
to be performed at temperatures of several /iK, an or- 
der of magnitude hotter than the current 100 to 300 nK 
experiments. At these temperatures the density needed 
for collapse is so high that inelastic two and three body 
processes seriously restrict the lifetime of the condensate 
| p6[ . Thus it may become necessary to work with softer 
traps (larger L), where the line of collapse lies at lower 
densities (see Fig. ||). 

The formalism used here to discuss the collapse of a gas 
with attractive interactions also describes domain forma- 
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FIG. 3. Maximum number of condensed particles as a func- 
tion of temperature for the same parameters as in Fig. ^. 
Present experiments are performed in the shaded region. The 
inset shows the maximum total density for the same temper- 
ature range. 



tion in spinor condensates, and gives a nice qualitative 
understanding of experiments at MIT [l8| in which op- 
tically trapped ^^Na is placed in a superposition of two 
spin states. Although all interactions in this system are 
repulsive, the two different spin states repel each other 
more strongly than they repel themselves, resulting in 
an effective attractive interaction. The collapse discussed 
earlier becomes, in this case, an instability towards phase 
separation and domain formation. The equilibrium do- 
main structure is described in JlOt . Here we focus on the 
formation of metastable domains. 

The ground state of sodium has hyperfine spin F = 
1. In the experiments the system is prepared so that 
only the states |1) = = l,mF = 1) and |0} = \F = 
l,mp = 0) are involved in the dynamics. The effective 
Hamiltonian is then 



H= d 



2m 



+ y(r)V'lV. + ^4V']V'jV'.- (10) 



where (i = 0,1) is the particle destruction operator for 
the state summation over repeated indices is assumed. 
The effective interactions, gij = inh^aij /m, are related 
to scattering amplitudes ap^o and ap=2, corresponding 
to scattering in the singlet (Fi + F2 — 0) and quintuplet 
(Fi -f ^2 = 2) channel by IfT 



Sa 



a = flu — floi — aio — ap=2, 
Oil — aoo = {ap=2 — aF=o)/3. 



(11a) 
(lib) 



a = 2.75nm and Sa = 0.19nm. We intro- 
duce g = inh^a/m and Sg = Anh^Sa/m. In the mean 



Numerically, 



field approximation, the interaction in Eq. (^^ becomes 
a function of n,„=o and n, the density of particles in the 
state |0) and the total density, respectively: 
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which shows expHcitly the effective attractive interac- 
tion. Initially the condensate is static withdensity n = 
lO^'^cm"'^, and all particles in state |1). A radio fre- 
quency pulse places half of the atoms in the |2) state 
without changing the density profile. Subsequently the 
two states phase separate and form domains from 10 to 
50 fiui thick. The trap plays no role here, so we can 
neglect V{r) in Eq. ( p^ ) and consider a uniform cloud. 

Linearizing the equations of motion with an equal den- 
sity of particles in each state, we find two branches of ex- 
citations corresponding to density and spin waves , 
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— -^~gn— + 0{6g), (13a) 
2m / m 

^\ ^5gn^ + 0{{5gf). (13b) 
2m / 4m 



Since 5g > Q we see the appearance of spin excitations 
with imaginary frequency. The mode with the largest 
imaginary frequency grows most rapidly, and the width 
of the domains formed should be comparab le to the wave- 
length A of this mode. By minimizing Eq. ( |l3b[ ) one finds 
A = ^Jl-K jn 6a — lO^m, in rough agreement with the ob- 
served domain size. 
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